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Abstract. We formulate the the problem of persistent diffusion in three dimensions 
from the perspective of the Frenet-Serret equations. In contrast to one and two 
dimensional systems, in three dimensions persistent diffusion is, in general, a third 
order process. In this paper we derive a Fokker-Planck equation for the process and we 
calculate its effective diffusion constant. We also provide expressions for the asymptotic 
average displacement of the walk, as well as explicit expressions for the Fourier-Laplace 
transform of the correlations between the tangent, normal and binormal vectors of the 
motion. 
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Persistent random walks have been proposed as models for a wealth of physical and 
biological systems [DlalEliaElE] . In the simplest scenario, this process is characterized 
by the tendency to continne moving in the direction of the last step. In certain 
circnmstances, this tendency can be thonght of a kind of inertial effect that keeps the 
particle moving in the direction of motion [7]. On the other hand, in the context of 
biology, the tendency may arise from the fact that animals move in the direction in 
which they are facing [8]. 

For one dimensional systems the probability distribntion fnnction for the position of 
the walker independent of the direction it is moving, satishes the telegrapher’s eqnation 
in the continnous limit [7]. In two dimensions, the probability distribntion fnnction 
for the position of the walker can be shown to satisfy a telegrapher’s eqnation only 
approximately, and little is known for three and higher dimensions [9l [10]. However, 
in all cases, the net effect of persistence is that it indnces a short term memory in the 
motion in the sense that the direction of the path is closely related to previous directions 
over a correlation time, but in the long run, the transport process remains diffusive. 

An interesting generalization of this process is the case in which the tendency to 
move is not in the same direction as the previous step, but rather, at a hxed deviation 
from the previous direction. In two dimensions, this modihcation induces a tendency 
to “loop” if, say, the direction of each step is narrowly centered at a hxed, non zero 
angle from the previous step. This modihcation was shown to give rise to nonmonotonic 
behavior of the walker’s ehective dihusion constant as a function of the width of the 
distribution of angle changes |T0]. For very small widths, the walker would perform 
many loops that brought it back close to it’s initial position and, thus, giving rise to 
small dihusion. If, on the other hand, the distribution is very broad, then the persistence 
ehects are quickly lost and the process is essentially a normal random walk. However, 
between both limits, large loop segments can become uncorrelated before returning near 
the starting point, giving rise to a large ehective dihusion coefficient |inj . 

In this work we address the three dimensional case, namely we study the transport 
properties of persistent dihusion with curvature and torsion. This three dimensional 
motion induces a tendency for the trajectories to spiral, a behavior that has been 
observed in many animals “... helical trajectories are nearly universal for organisms 
that are less than 0.5 mm long and live in water.” [T] 

Further, in spite of its limitations, persistent random walks are useful models for 
the shape of linear polymers with some degree of bond rigidity; thus persistent dihusion 
with curvature and torsion would be appropriate models for polymers that tend to form 
helices jH [6] . 

As in [5], we formulate the problem in terms of the Frenet-Serret equations with 
noise terms added to the curvature and torsion. We then construct a Fokker-Planck 
equation m for this process, from which we can calculate the correlation functions of 
the Frenet frame. In particular, we can evaluate the velocity autocorrelation function 
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from which the effective diffusion constant of the process is obtained. 
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2. Formalism 


To formulate the persistent diffusive process, we consider the Frenet description of a 
trajectory [DEI El [5], which consists in the evolution of the right handed frame of unit 
vectors given by T, the vector tangent to the curve, the normal N, in the direction of 
the derivative of T, and the binormal B = T x N. The Frenet-Serret equations then 
state that 
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where s is the arclength, K and T are the curvature and torsion of the curve respectively, 
which can, in principle, depend on s too. The time evolution of the particle following 
the trajectory is given by 

dX 
dt 

where v is the speed of the particle, which we will take as a constant throughout this 
work.. 

To simplify the construction of the Fokker-Planck equation, we begin by discretizing 
the Frenet-Serret equation as follows. We write the arclength as s = nXs and take the 
following discretized version of the evolution equation: 
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The rationale for this choice is that at each step, the Frenet frame undergoes a pure 
rotation. Equivalently, inverting eq.jjS]) we have 
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Finally, we introduce noise by taking 

Kn —>■ kXs + air]i{n){XsY^‘^ 
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Where k and r will be assumed to be constant, while rii{n) and r] 2 {n) are random 
variables satisfying 

{r]j{n)) = 0 {r]j{n)r]i{m)) = Sij6n,m (7) 

To order As we can evaluate the matrix as follows 
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where the n-dependencies have been suppressed from the last expression in the interest 
of typographical clarity. 

Thus, if we denote Pn{T, N, B) the probability density function for the Frenet frame 
being given by the vectors T, N and B at step n, then, correct to order As, we derive 
the Fokker-Planck equation from the relation 

Pn+Yf, N, B) = l^Pn , A, B)) ) . (9) 

Here the angular brackets denote the average over the random noise r]i{n) and ri 2 {n) 
and R~^ is as in (0). 

To lighten the expressions to come, let us now dehne the drifts 
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Putting in ([9]), PniT, N, B) = P(T, A, B] s), with s = nAs, expanding to linear order 
in As and hnally dividing by As, we obtain 
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where Vt, V n and V b denote the gradient in the T, A and B coordinates respectively. 
The equation above describes the evolution of the probability distribution for the 
components of the Frenet frame along the arclength of the particle’s trajectory and 
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becomes exact in the limit As —>■ 0. Equation fllip can be rewritten in a slightly more 
symmetric form as 
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either way, the equation seems to be rather intractable. However, it can be used to 
study the evolution of the moments of T, N and B] from these, the transport properties 
of x{t) can be inferred. In particular, by carrying out the integrations, one hnds that 
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Here the dependence on the initial condition has been written explicitly to emphasize the 
fact that these statistics give us access to the time correlation functions. The solution 
can be written as 
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where T is the matrix appearing in eq. flT^ . Thus, the matrix of correlation functions 
can be written explicitly as 

/ it ■ To) it ■ No) {t ■ Bo) \ 

M(s)= (4-To) (4-No) {4-Bo) (15) 

\ {B, ■ To) (4 • No) (4 • Bo) 

The Laplace-Fourier transform of the correlation matrix is thus given by 

M{u) = {iu - T)-\ (16) 

If we now dehne A(a;) as follows 

A(a;) = (4+2*C(;)((j2+2^C(;)(4+4+2*<^)+4A;^((j2+2*tu)+4r^(cr42fC(;), (17) 

we find the matrix A{u)Ai{u))/2 to be given by 
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Note that the matrix is not symmetric as the process is not invariant under simple 
time inversion. Rather, T and B are odd, whereas N is even, leading to the observed 
alternation between symmetric and antisymmetric elements in the matrix. 
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The explicit inversion of the above results to obtain the time-dependent correlation 
functions can, of course, be carried out in principle, as it only requires the solution 
of a third order polynomial. Nevertheless, the results are messy and unenlightening. 
Instead, we focus on particular cases which illustrate different aspects of the process. 
One simple situation is reached in the planar case, in which both the torsion r and its 
associated noise vanish (<72 = 0). In this case, the binormal vector is constant and the 
motion is restricted to a plane. In this situation, ((18 
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where A(a;) here reduces to 

A2d(n;) = 2iuj[{al + 2iujf + Ak\ (20) 

The s-dependent correlations can also be determined explicitly in this case. Indeed one 
readily hnds 
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which corresponds to circular motion which is damped due to the noise. 

Further, we can use AA to write the Fourier-Laplace transforms for the average of 
the Frenet vectors; for example, for the mean tangent vector we have: 
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This quantity connects to the behavior of the walker’s average position by 
{xi^\TotNo,Bo) = —fAuj\To^No,Bo). 

lU 

Thus, in the long time limit, the average position of the persistent walker saturates at 
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Another interesting quantity is the Darboux vector [T], which in our case, will be writen 
as A = tT A- kB. This vector is akin to the angular momentum of the walker. In the 
absence of noise, this vector is constant, as can be seen directly from flT5]l . and lies on 
the axis of the helical motion. In the presence of noise, however, the Darboux vector 
evolves according to 
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Here = tTq + kB^ is the initial value of A. It is worth noting that in the small noise 
limit, the saturation value of the average position, given in fl2^ is given by 

2tAo 
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(26) 


which reflects the intuitively appealing fact that in the low-noise limit, the net 
displacement is in the direction of the axis of the helix. 

Perhaps more interestingly, from the correlation functions we can compute the 
effective diffusion constant. Assuming, again, that the instantaneous speed v is constant, 
we can write 
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recalling that s* = vti. The effective diffusion constant D will be given by 


D = lim 

t^OO 


{x{t)-x{ t)) 
2t 


= V 


To • Tn ) ds = 2v 


+ cr|) + 4 t’ 


cr^(7|(cr^ + <^l) + + 4r^al 


(27) 


(28) 


This expression’s behaviour becomes easier to visualize after introducing scaled versions 
of the curvature, the torsion and the diffusion constant as well as the adimensional 
parameter k, which measures the relative importance of the curvature and the torsion: 
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Thus D grows with r and decreases with k. Intuitively, this corresponds to the fact 
that torsion favors a rectilinear motion in the direction of the Darboux vector, thus 
favouring diffusion by enhancing transport, whereas curvature leads to looping around 
in circles, which conhnes the walker’s motion. For the true diffusion constant, matters 
are a bit more complex, since it depends really on three parameters. Thus, for instance, 
upon increasing Ui at hxed af -\- —that is, at hxed total noise—we decrease k and 

thus increase D. The behaviour of D, however, is in general non-monotonic, since it 
corresponds to D divided by a^. 

It is interesting to note that fl28|) recovers the two-dimensional result taking the 
torsion r —)■ 0 hrst, simplifying and then <72 —)■ 0. This procedure yields 
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which is analogous to the result obtained in m. We note that, as above, D 2 d is a non 
monotonic function of Ui. It is also interesting to note that a different, intrinsically 
three dimensional result, is obtained if we take the limits in opposite order. Indeed, if 
we let ct 2 —>■ 0 first, we see that k —>■ 0 and therefore 

Dh = —V (32) 

which is independent of all the other parameters of the problem. From the scaling 
form fl30|) . it follows quite straightforwardly that the same result also holds if /c = 0. 
Intuitively, however, it is quite far from being obvious why, in the former case, the 
diffusion constant only depends on ai and v. 

Indeed, it is quite clear that a non-trivial crossover behaviour must be observed 
when we have a 2 = 0 and r small. In this case, we might have expected the behavior 
to be close to that for r = 0, yet the exact expressions for the diffusion constants differ 
by a finite amount. This means that, in this case, we shall observe two-dimensional 
diffusion for a long-time, and would expect to see the corresponding diffusion constant, 
see fl3T]l over a limited time range. After that, we would expect torsion to bring the 
system significantly out of the plane, and thus to lead to the expression fl52]) . Indeed, 
this can be verified by a straightforward though tedious analysis of the zeroes of A(a;). 
One finds a low-lying imaginary zero with the value 

2ia‘lT‘^ 

^ af + 4fc2 

which corroborates the picture 

When A; = 0, on the other hand, the tendency of the process is to move, on average, 
in the same direction of the previous motion, so that we have a typical persistent walk. 
The torsion, on the other hand, with or without noise, merely rotates the new step 
around the previous tangent and thus has no influence on the diffusive properties of the 
process. 



4. Asymptotic shape of the average trajectory in the low-noise limit 


As a final remark, we observe a peculiar transition in the large-time behaviour of the 
average trajectory. Indeed, the average position of the trajectory at arclength s is 
determined by the matrix exp(rs). For large s it must therefore be dominated by that 
eigenvalue of F having the smalles real part in absolute value. As is readily seen, the 
three eigenvalues of F in the limit of small noise are given by 
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up to terms of higher order in the noise. Thus, in the low-noise limit, we have two 
groups of eigenvalues, one tending to 0 and corresponding essentially to the Darboux 
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vector as an eigenvector. This thus describes the motion in the direction of the helix 
axis. The other two eigenvalues correspond to eigenvectors determining the rotational 
motion. As we see, their real parts are not equal, and an easy calculation shows that 
Ao dominates if 

$ = (2A;2 - T^)al + (2 t^ - > 0 (35) 

and A-t otherwise. Thus, according to the sign of <h, we have two qualitatively different 
motions: if $ > 0, the circular motion decays more rapidly than the progressive motion, 
and the asymptotic shape of the particle’s trajectory is a very loosely wound helical 
motion, in which the particle in one rotation advances by a distance far larger than 
the helix’ radius. If < 0, on the other hand, the opposite situation prevails, and the 
motion is on a circle of a radius much larger than the distance by which the particle 
advances for each rotation. In other words, we might say that the asymptotic pitch of 
the average motion diverges at large times for $ > 0, whereas it vanishes for $ < 0. We 
show in Figure 1 the autocorrelation function of the tangent vector T for large times 
for 2 cases in which <F > 0 and <F < 0 respectively. The difference is quite striking. 

Finally, let us point out that the above-mentioned transition does not always occur. 
In particular, if k and r are equal, or of similar size, then we always have <F > 0. For 
both signs of <F to be possible, we require 

max ^ > \/2. (36) 

When $ = 0, the model maintains a truly helicoidal dynamics for all times. It would 
be of interest to know whether any particular significance attaches to this fact in real 
systems . 
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Fig. 1: The autocorrelatin function for T for k = 0.4 and r = 1 in the cases where 
(Ti = (72 = 0.1 (upper hgure) and cxi = 0.15, a 2 = 0.05 (lower hgure). 


5. Conclusions 

Let us briefly summarize what we have done: we have dehned a discretization of the 
Frenet-Serret equations describing a curve in three-dimensional space. Assuming that 
the curvature and torsion are both constant, and proportional to the discretization 
parameter As this equation we have dehnes a discretized helical motion. We add 
uncorrelated Gaussian noise of the order of (As)^/^ to the curvature and torsion to dehne 
a helical brownian motion. Next derive the Fokker-Planck equation for this stochastic 
process. The resulting process is, in fact, a random walk in which the probability of 
every step depends on the step’s relative position with respect to the preceding two 
steps. Indeed, the curvature dictates the average angle to the preceding step, whereas 
the torsion determines the average angle with respect to the plane generated by the two 
earlier steps. The process is fully determined by hve parameters: the average curvature 
A:, the average torsion r, the noise intensity for the curvature af, the noise intensity 
for the torsion cr| and hnally the speed v at which the particle moves. While this last 
parameter only hxes an otherwise arbitrary time scale and can thus in principle be 
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disregarded, the other parameters are all essential. 

We then proceeded to determine the Fourier-Laplace transform of the various 
correlation functions of the orthogonal triad T, N and B. From this it is possible to 
determine the effective diffusion constant. An explicit formula for it is derived in terms 
of the model’s hve parameters, see fj28l) . Interestingly, there is a scaling parameter k 
dehned in terms of these hve, see ([3]), such that an appropriately rescaled diffusion 
constant depends only on k, see fl30|) . 

Other quantities which can be evaluated explicitly are the asymptotic average 
position (as a vector) of the walk, as well as the correlation function of the Darboux 
vector, which corresponds to the axis of the helix which would be generated in the 
noiseless limit. We show in particular that the asymptotic average position of the walk 
indeed coincides with the Darboux vector in the low-noise limit. 

Various special cases can also be treated more or less explicitly. Among these are 
the two-dimensional case ((T 2 = 0 and r = 0), for which we have given an explicit 
expression for the correlations. 

A still open question concerns the response of the process when an external held is 
applied; the difficulties for advancing in this question arise from the fact that a spatially 
hxed held rotates in the Frenet frame. A related problem would be to understand 
the ehects in the transport properties due to periodic or strategic variations in the 
parameters, attempting to steer the particle in a given direction [DEI El [5]. 
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